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Abstract. After a brief discussion on neutrino electromagnetic properties, we consider the 
problem of neutrino energy spectra in different media. It is shown that in two particular cases 
(i.e., neutrino propagation in a) transversally moving with increasing speed medium and b) 
rotating medium) neutrino energies are quantized. These phenomena can be important for 
astrophysical applications, for instance, for physics of rotating neutron stars. 



1. Introduction - neutrino electromagnetic properties 

Initially the problem considered in this paper has originated from the studies of neutrino 
electromagnetic properties and related items. There is no doubt that the recent experimental and 
theoretical studies of flavour conversion in solar, atmospheric, reactor and accelerator neutrino 
fluxes give strong evidence of non-zero neutrino mass. A massive neutrino can have non-trivial 
electromagnetic properties [1]. A recent review on neutrino electromagnetic properties can be 
found in [2]. The present situation in the domain is characterized by the fact that in spite of 
reasonable efforts in studies of a neutrino electromagnetic properties, there is no any experimental 
confirmation in favour of a neutrino electromagnetic characteristics being nonvanish. However, 
it is very plausible to assume that a neutrino may have nonzero electromagnetic properties. In 
particular, it seem very reasonable that neutrino have nonvanishing magnetic moment [2,3]. 

Neutrino magnetic moment interaction effects. If a neutrino has non-trivial electromagnetic 
properties, notably non- vanishing magnetic (and also electric (transition) dipole moments or non- 
zero millicharge and charge radius), then a direct neutrino couplings to photos become possible 
and several important for applications processes exist [4]. A set of typical and most important 
neutrino electromagnetic processes involving the direct neutrino couplings with photons is: 1) 
a neutrino radiative decay z^i — > i/2 + 7, neutrino Cherenkov radiation in external environment 
(plasma and/or electromagnetic fields), 2) photon (plasmon) decay to a neutrino-antineutrino pair 
in plasma 7 vv, 3) neutrino scattering off electrons (or nuclei), 4) neutrino spin (spin-flavor) 
precession in magnetic field. Another very important phenomena is the resonant amplification of 
the neutrino spin- flavour oscillations in matter that was first considered in [5]. 
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Note a new mechanism of electromagnetic radiation produced by a neutrino moving in matter 
and originated due to neutrino magnetic moment [6-8]. It was termed the spin hght of neutrino 
in matter (SLv) [6]. Although the SLv was considered first within quasiclassical approach, it was 
clear that this is a quantum phenomenon by its nature. The quantum theory of this radiation 
has been elaborated [7] (see also [8]) within development [9, 10] of a quite powerful method that 
implies the use of the exact solutions of the modified Dirac equation for the neutrino wave function 
in matter. For elaboration of the quantum theory of the SLv one has to find the solution of the 
quantum equation for the neutrino wave function and for the neutrino energy spectrum in medium. 



2. Quantum equation for neutrino in medium 

The modified Dirac equation for the neutrino wave function exactly accounting for the 
neutrino interaction with matter [7]: 

{n^.^'^-l7^.i^ + l5)F-rn]^i^)^0■ (1) 

This is the most general form of the equation for the neutrino wave function in which the 
effective potential = ^(1 + 75)//^ includes both the neutral and charged current interactions of 
neutrino with the background particles and which can also account for effects of matter motion 
and polarization. It should be mentioned that other modifications of the Dirac equation were 
previously used in [11] for studies of the neutrino dispersion relations, neutrino mass generation 
and neutrino oscillations in the presence of matter. 

In the case of matter composed of electrons, neutrons, and protons and for the neutrino 
interaction with background particles given by the standard model supplied with the singlet right- 
handed neutrino one has 

r^V2GF fhf + ^'<if\ (2) 



f=e,p,n 



where 



,W = (4{)-2Q(/)sin^^^^+<5.,),,f =-(4{'+M,'^e/ = { \J;ii^;_ (3) 



Here 1^^, ^^^'^ Q^'^^ ^'^'3, respectively, values of the isospin third components and the electric charges 
of matter particles (/ = e,n,p). The corresponding currents j'^ and polarization vectors are 

j; = («/,«/v/), - f»/(C/v/),nK/yr^+ ''f^'f'^^IllA , (4) 

where 0w is the Weinberg angle. In the above formulas (g]), Uf, v/ and Cf (0 < IC/I^ < 1) stand, 
respectively, for the invariant number densities, average speeds and polarization vectors of the 
matter components. 

In the case of matter at rest it is possible to solve the modified Dirac equation for different 
types of neutrinos moving in matter of different composition, as it is shown in [7]. The energy 
spectrum of different neutrinos moving in matter is given by 



= erj^p^^l — sa—"^ + + am. (5) 

In the general case of matter composed of electrons, neutrons and protons the matter density 
parameter a for different neutrino species is 

1 Gf { \ 
= -^-^ — (^neiisin^ 9w + g) + np{l - 4sin^6li4') - n„j, (6) 

where g = I for the electron neutrino and g = —1 for the muon and tau neutrinos. 

The value rj =sign(l — say) in ^ provides a proper behavior of the wave function in the 
hypothetical massless case. The values s = ±1 specify the two neutrino helicity states, p+ and v-. 
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The quantity e — ±1 splits the solutions into the two branches that in the limit of the vanishing 
matter density, a ^ 0, reproduce the positive- and negative-frequency solutions, respectively. 

In the next two sections we apply the developed method of exact solutions to two particular 
cases when neutrino is propagating in transversally moving with increasing speed medium [9, 12] 
and a rotating medium of constant density. In both cases the obtained energy spectrum of neutrino 
is quantized like the energy spectrum of an electron is quantized in a constant magnetic field. 



3. Neutrino quantum states in transversally moving with increasing speed medium 

First we consider a neutrino propagating in medium composed of neutrons that move 
perpendicular to the neutrino path with linearly increasing speed. This can be regarded as the 
first approach to modelling of neutrino propagation inside a rotating neutron star [9,12]. The 
corresponding modified Dirac equation for the neutrino wave function is given by ([T]) with the 
matter potential accounting for rotation, 

r = -G(n,nv), v=(c^y,0,0), (7) 

where G = Here uj is the angular frequency of matter rotation around OZ axis, it also is 

supposed that the neutrino propagates along OY axis. For the neutrino wave function components 
i^{x) we get from the modified Dirac equation ([1]), a set of equationfQ, 



(^0 - 93) + Gn] V^i + [ - {idi + 32) + Gnujy 
[ i^tdi + 82) + GnLuy]^^ + [i {do + 83) + Gn 



-02 = mip3, 
-02 = rmpi, 



(8) 

i (do + 83) -03 + {idi + 82)^4 = rmpi , ^ 
{i8i - 82) -03 + i{8o- 83) -04 = rmp2. 

In general case, it is not a trivial task to find solutions of this set of equations. 

The problem is reasonably simplified in the limit of a very small neutrino mass, i.e. when the 
neutrino mass can be ignored in the left-hand side of ([H]) in respect to the kinetic and interaction 
terms in the right-hand sides of these equations. In this case two pairs of the neutrino wave 
function components decouple one from each other and four equations ([H]) split into the two 
independent sets of two equations, that couple together the neutrino wave function components 
in pairs, (tAi, -02) and (V'3, ^4)- 

The second pair of equations ([5]) does not contain a matter term and is attributed to the 
sterile right-handed chiral neutrino state, -0^. The corresponding solution can be taken in the 
plain-wave form 

0_R ~ exp{i{~pot + pix + p2y + Psz)}^^, (9) 

where is the neutrino momentum. Then for the components -03 and -04 we obtain from ([H]) the 
following equations 

{Po - Ps) -03 - {Pi - iP2) = 0, 
- {pi + ip2) -03 + (po + Pz) ^^4 = 0. 

Finally, from (flO)) for the sterile right-handed neutrino we get 



0fl 



/ \ 





L3/2^2po(po-P3) 



(11) 



~Pl + W2 

V Pz-Po J 

where px — p^^x'', p^ = {pa,Pi,P2,P3) and — {t, x, y, z). This solution, as it should do, has the 
vacuum dispersion relation. 

In the neutrino mass vanishing limit the first pair of equations ([8]) corresponds to the active 
left-handed neutrino. The form of these equations is similar to the correspondent equations for a 
charged particle (e.g., an electron) moving in a constant magnetic field B given by the potential 
A — {By, 0, 0) (see, for instance, [13]). To display the analogy, we note that in our case the matter 



X The chiral representation for Dirac matrixes is used. 
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current component nv plays the role of the vector potential A. The existed analogy between an 
electron dynamics in an external electromagnetic field and a neutrino dynamics in background 
matter is discussed in [9]. 

The solution of the first pair of equations ([5]) can be taken in the form 

ipL'^T exp{i(-poi + Pix + P3z)}ip{y), (12) 



and for the components -01 and -02 of the neutrino wave function we obtain from 
equations 



0, 



the following 



(13) 



where 



p = Gnu. 



For the wave function we finally get 

-ipot+ipix+ip^z 



( iPO 



■P3 



L^(po-p3 + Gny + 2pN 



Gn)uN{j]) \ 



V 



-V2piVujv_i(7?) 





(14) 



(15) 



where U]\[{r]) are Hermite functions of order N. For the energy of the active left-handed neutrino 
we get 



Po 



pl + 2pN - Gn, TV = 0,1,2,.. 



(16) 



The energy depends on the neutrino momentum component p^ along the rotation axis of matter 
and the quantum number N that determines the magnitude of the neutrino momentum in the 
orthogonal plane. For description of antineutrinos one has to consider the "negative sign" energy 
eigeinvalues (see similar discussion in Section 2.2). Thus, the energy of an electron antineutrino 
in the rotating matter composed of neutrons is given by 



Po 



pl + 2pN + Gn, TV = 0,1,2,, 



(17) 



Obviously, generalization for different other neutrino flavours and matter composition is just 
straightforward (see ([4]) and 

Thus, it is shown [12] that the transversal motion of an active neutrino and antineutrino is 
quantized in moving matter very much like an electron energy is quantized in a constant magnetic 
field that corresponds to the relativistic form of the Landau energy levels (see [13]). 



4. Neutrino energy in a rotating medium 



Now we consider a more consistent model of a neutrino motion in a rotating matter. For this 
case we choose the effective matter potential in ([T]) in the following form 

= -G(n, nv), V = (-wy, ux, 0). (18) 

Contrary to the case considered in the previous section, equation ([1]) with the potential (jlSp 
describes the case when a neutrino is moving in a rotating media. It is shown below how the 
equation ([1]) with (|18p can be solved and the corresponding neutrino energy spectrum is obtained. 
The solution of the equation ([T]) with (fT8|) can be sought in the form 

/ V'i(a;,y) \ 

^3{x,y) 
\ i^i{x,y) ) 



il;{t,x,y,z) = ^ 



(19) 
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Substituting (|19p into JT]) with (fT5|) and using the expHcit form of the Dirac matrices in the 
chiral representation, we arrive at a system of hnear equations for the neutrino wave function 
components: 



(20) 



(21) 



-(po +P3 + Gn)^i + i - + Gnuj{x - iy)| -02 = -"ii/'s 

* { ( + ~ Gnuj{x + ij/)| -01 + (j33 - Po - Gn)?/;2 = -m-0, 

(po - P3)'03 + « - ^"4 = TO0 

+* + ^3 + (PO + P3)'04 = ™0: 

In the polar coordinates x + iy = re*"^, x — iy = re"**^ one has 

_9_ I — ( -I- 1^ ) M. — e~*"^ ( — i-^ 1 

Ox c?y \ dr r dtj) J ' dx dy \ dr r d(p J ' 

and the system of equations ([^0)1 transforms to 

-(po +P3 + Gn)^i + ze-'* - + pr} 02 = -"^03, 

(po -P3)V'3 + *e-** - 1^) ^4 = r7i0i, 

+^e'* + 7w) ^3 + (Po +P3)V'4 = m^2. 

It is possible to show that the operator of the total momentum Jz — Lz+Sz, where Lz = ~i-§^^ 
Sz = \o'3j commutes with the corresponding Hamiltonian of the considered system. Therefore the 
solutions can be taken in the form 



(23) 



that are the eigenvectors for the total momentum operator Jz with the corresponding eigenvalues 
Z — i. After substitution of the system ([2^ can be rewritten in the following form, 



(22) 



( ^ 
















iX3{r)e'^^~^^^ 


V ^4 ) 




\ X4(r)e^'* J 



(24) 



-(Po +P3 + Gn)xi + {-^ + i + Pr} X2 = -"^X3, 

{5^ ' ^ ^ ^'^l + (Po - P3 + Gn)x2 = mx4, 

(Po -P3)X3 + + ^) X4 = mxi, 

(-^ ~ ^) X3 - (Po +P3)X4 = -mX2- 

For the further consideration it is convenient to introduce the rising and decreasing operators 

d I — 1 „ d I , ^ 

pr, R = — H V pr. (25) 



(26) 



dr r dr r 

After applications of the decreasing R~ and increasing i?+ operators to the second and first 
equations of ([22]) correspondently, one gets 

R^R+Xi + ((Po + Gnf -pI- m^^ xi = m{Gnx3 + prxi), 
R^R^X2 + ((Po + Gnf -pI~ m^) X2 = m{Gnx^ + prX3)- 

Note that the system ([24|) . as well as the system (|8]), can be solved exactly in the limit of 
vanishing neutrino mass to — s- 0. In order to find the non-zero-mass correction to the energy 
spectrum of a neutrino in a bound state in matter, the neutrino square integrable wavefunction 
should be found. Therefore, in analogy with the zero-mass case we take X3 = X4 = in the lowest 
order of perturbation series expansion. Thus we arrive to the system 



i?-i?+Xi + ((po + Gn)2 
i?+i?"X2+ (po + Gn)2 



pi - Xi 

Pi 



0, 

TO^) X2 = 0. 



(27) 
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The solution of (I27p can be written in the form 
Xi\_( C,Ci-Hpr') 



X2 I \ C^C'sipr') 



(28) 



where C\ are the Laguerre functions [13]. After substitution of to ([?/)) and taking into account 
the properties of the increasing and decreasing rising operators, 

i?+4-i(pr2) = -2^^(7+0 4(pr2), 

R-C\{pr^) = 2VK^4-Hpr-2), ^ ^ 

we get from (j27p the equation for the neutrino energy spectrum in matter 

+pI+A{s + 1)p- [pa + Gnf =Q. (30) 

Solving this equation we get for the neutrino energies 



Po = ±Vm2+p2 + 4iVp-Gn, iV = 0, 1, 2, ... . (31) 

where the quantum number N = s + / is introduced. As usually, two signs in the solution 
corresponds to the neutrino and antineutrino energies, correspondently, 



Pq = ^m? +pI + ANp - Gn, pQ = ^m? +pl + AN p + Gn, . (32) 

From the obtained energy spectrum it just straightforward that the transversal motion 
momentum of an antineutrino is given by 



pi_ = 2VNGuj. (33) 

The quantum number N determines also the radius of the quasiclassical orbit in matter (it is 
supposed that N ::§> 1 and ps = 0), 



It follows that antineutrinos can have bound orbits inside a rotating star. To make an estimation 
of magnitudes, let us consider a model of a rotating neutron star with radius Rns = 10 km, matter 
density n = lO^^cm"'^ and angular frequency u — 2tt x 10^ s^^. For this set of parameters, the 
radius of an orbit is less than the typical star radius Rns if the quantum number TV < N^ax = 10^*^. 
Therefore, antineutrinos that occupy orbits with TV < 10^" can be bounded inside the star. The 
scale of the bounded antineutrinos energy estimated by ([321) is of the order po ^ ^ eV. It should be 
underlined that within the quasiclassical approach the neutrino binding on circular orbits is due 
to an effective force that is orthogonal to the particle speed. Note that there is another mechanism 
of neutrinos binding inside a neutron star when the effect is produced by a gradient of the matter 
density (see the last paper in [11]). A discussion on the "matter- induced Lorentz force" that can 
be introduced in order to explain a neutrino motion on quasiclassical circular orbits can be found 
in [9]. 

We ague that the effect of a neutrino energy quantization can have important consequences 
for physics of rotating neutron stars. 
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